%obj = KPOINTS('K-points', 0, 'Monkhorst Pack', [8 8 8; 0 0 0]);

clear;
latticeConstant = 5.474; % Angstroms
alpha = 55.195; % degrees


% dummy = [latticeConstant*cos(alpha*degree/2) -latticeConstant*sin(alpha*degree/2) 0;...
%   latticeConstant*cos(alpha*degree/2) latticeConstant*sin(alpha*degree/2) 0;...
%  latticeConstant*cos(alpha*degree)/cos(alpha*degree/2) 0 latticeConstant*sqrt(1-cos(alpha*degree)^2/cos(alpha*degree/2)^2)];
% 


lo = Lattice.create([latticeConstant latticeConstant latticeConstant], [alpha alpha alpha]*degree);
lo.transform_standard_rhombohedral;


obj = KPOINTS.create_RHL_1_bandstructure(alpha*degree);
obj.write;